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Abstract 

We present the dominant two-loop electroweak corrections to the partial decay 
widths to gluon jets and prompt photons of the neutral CP-odd Higgs boson A , 
with mass M^o < 2Myp, in the two-Higgs-doublet model for low to intermediate 
values of the ratio tan/3 = i^/vi of the vacuum expectation values. They apply as 
they stand to the production cross sections in hadronic and two-photon collisions, 
at the Tevatron, the LHC, and a future photon collider. The appearance of three 
75 matrices in closed fermion loops requires special care in the dimensional regu- 
larization of ultraviolet divergences. The corrections are negative and amount to 
several percent, so that they fully compensate or partly screen the enhancement due 
to QCD corrections. 
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The search for Higgs bosons is among the prime tasks at the Fermilab Tevatron and will 
be so at the CERN Large Hadron Collider (LHC), to go into operation later during this 
year, and the International e + e~ Linear Collider (ILC), which is currently being designed. 
The standard model (SM) contains one complex Higgs doublet, from which one neutral 
CP-even Higgs boson (H) emerges in the physical particle spectrum after the electroweak 
symmetry breaking. Despite its enormous success in describing almost all experimental 
particle physics data available today, the SM is widely believed to be an effective field 
theory, valid only at presently accessible energy scales, mainly because of the naturalness 
problem related to the fine-tuning of the cut-off scale appearing quadratically in the Higgs- 
boson mass counterterm, the failure of gauge coupling unification, the absence of a concept 
to incorporate gravity, and the lack of a cold-dark- matter candidate. Supersymmetry 
(SUSY), which postulates the existence of a partner, with spin shifted by half a unit, to 
each of the established matter and exchange particles, is commonly viewed as the most 
attractive extension of the SM solving all these problems. The Higgs sector of the minimal 
SUSY extension of the SM (MSSM) consists of a two-Higgs-doublet model (2HDM) and 
accommodates five physical Higgs bosons: the neutral CP-even h° and H° bosons, the 
neutral CP-odd A boson, and the charged if ± -boson pair. At the tree level, the MSSM 
Higgs sector has two free parameters, which are usually taken to be the mass M^o of the 
A boson and the ratio tan/3 = V2/V1 of the vacuum expectation values of the two Higgs 
doublets. 

The discovery of the A boson would rule out the SM and, at the same time, give 
strong support to the MSSM. At the LHC, this will be feasible except in the wedge of 
parameter space with M^o ;> 250 GeV and moderate value of tan /3, where only the h° 
boson can be detected pQ. For low to intermediate values of tan/3, gluon fusion is by 
far the dominant hadroproduction mechanism. At large values of tan f3, A°bb associated 
production becomes important, too, especially at LHC cm. energy, -^i = 14 TeV [2J. At 
the ILC operated in the 77 mode, via Compton back-scattering of highly energetic laser 
light off the lepton beams, single production of the A boson will allow for its discovery, 
also throughout a large fraction of the LHC wedge, and for a precision determination of 
its profile [3J. Two-photon collisions, albeit with less luminosity, will also take place in the 
regular e + e~ mode of the ILC through electromagnetic bremsstrahlung or beamstrahlung 
off the lepton beams. 

In the mass range M40 < 2m t and for large values of tan f3 in the whole M^o range, 
the A boson dominantly decays to a bb pair, with a branching fraction of about 90% [2]. 
As in the case of the H boson of the SM, the rare 77 decay channel may then provide 
a useful signature at the LHC if the b and b quarks cannot be separated sufficiently well 
from the overwhelming background from quantum chromodynamics (QCD). The A — > gg 
channel will greatly contribute to the decay mode to a light-hadron dijet, which will be 
measurable at the ILC. 

Since the A boson is neutral and colorless, the t4°77 and A°gg couplings are loop 
induced. As the A boson has no tree-level coupling to the W boson and its coupling 
to sfermions flips their "handedness" (left or right), the v4°77 coupling is mediated at 
leading order (LO) by heavy quarks and charged leptons and by light charginos [I]. The 
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A°gg coupling is generated at LO by heavy-quark loops [5]. 

Reliable theoretical predictions for the A 0, yy and A°gg couplings, including higher- 
order radiative corrections, are urgently required to match the high precision to be reached 
by the LHC and ILC experiments [6]. Specifically, the properties of the A boson, es- 
pecially its CP-odd nature, must be established, and the sensitivity to novel high-mass 
particles circulating in the loops must be optimized. The present state of the art is 
as follows. The next-to-leading-order (NLO) QCD corrections, of relative order 0(a s ) 
in the strong-coupling constant a s , to the partial decay widths T(A° — > 77) E] and 
T(A° — > gg) [8], and the production cross section a(gg — > A ) [9j E] are available for 
arbitrary values of quark and A°-boson masses as one- dimensional integrals, which were 
solved in terms of harmonic polylogarithms for T(A° — > 77), T(A° — > gg), and the virtual 
correction to cr(gg — > A ) [Till ITT] . The latter was also obtained for general color factors 
of the gauge group SU(iV c ) in the limit m t — > 00 using an effective Lagrangian [12] . The 
next-to-next-to-leading-order (NNLO) QCD corrections, of 0(al), to T(A° — > gg) [T3] 
and <r(gg — > A ) [FT] were found for m t — > 00 using an effective Lagrangian. The 0(a s ) 
SUSY QCD correction, due to virtual squarks and gluinos besides the heavy quarks, to 
cr {99 ~~ * ^4°) was obtained from an effective Lagrangian constructed by also integrating out 
the SUSY particles [15]. The two-loop master integrals appearing in the latter calculation 
if the masses of the virtual scalar bosons and fermions are kept finite were expressed in 
terms of harmonic polylogarithms [TT] . 

In this Letter, we take the next step and present the dominant electroweak corrections 
to T(A° — > 77) and T(A° — > gg) at NLO. Since these are purely virtual, arising from 
two-loop diagrams, they carry over to (7(77 — > A ) and u{gg — > A ), via 

°{rtl99 - A ) = ^^r(A° -> vy/gg)* (s - M%) , (1) 

where iV 7 = 1 and N g = N% — 1 = 8 are the color multiplicities of the photon and the 
gluon, respectively, and s is the partonic cm. energy square. For the time being, we 
focus our attention on the particularly interesting region of parameter space with low to 
intermediate Higgs-boson masses, M h o, M H o, M A o, M H ± < m t , and low to moderate value 
of tan j3, tan (3 <C rrit/mb, and assume that the SUSY particles are so heavy that they can 
be regarded as decoupled, yielding subdominant contributions. The dominant electroweak 
two-loop corrections are then of relative order O(xt), where Xt = G fvo^ j (8ti 2 \/2) w 
3.17 x 10~ 3 with Gf being Fermi's constant. In the case of T(A° — > 77), they arise from 
the class of generic Feynman diagrams shown in Fig. HI which, besides the t and b quarks, 
involve the charged and neutral Goldstone bosons, and z°, and the five Higgs bosons. 
Here, it is understood that the b quark only couples to the w ± and H ± bosons because its 
couplings to the neutral scalar bosons are suppressed unless tan/3 <^m t /m b and, of course, 
that the z°, h°, H°, and A bosons do not couple to the photon. We explicitly checked 
that the and Z° bosons do not contribute at 0(xt). The diagrams contributing to 
T(A° —>■ gg) at 0(x t ) emerge from Fig. CD by omitting those where a photon couples to a 
scalar boson and by replacing the photons by gluons. 
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Figure 1: Feynman diagrams contributing to T(A° — ► 77) at 0(xt). S = 
w^, z°, h°, H°, A , if ± and f = t,b denote generic scalar bosons and fermions, respec- 
tively. The couplings of the z°, h°, H°, and A bosons to the b quark are to be neglected 
and those to the photon vanish. 

We now outline the course of our calculation and exhibit the structure of our results. 
Full details will be presented in a forthcoming communication [16]. Since we consider the 
SUSY partners to be decoupled, we may as well work in the 2HDM without SUSY. We 
may thus extract the ultraviolet (UV) divergences by means of dimensional regularization, 
with D = 4 — 2e space-time dimensions and 't Hooft mass scale \i. For convenience, 
we work in 't Hooft-Feynman gauge. We take the Cabibbo-Kobayashi-Maskawa quark 
mixing matrix to be unity, which is well justified because the third quark generation is, 
to good approximation, decoupled from the first two. We adopt Sirlin's formulation of 
the electroweak on-shell renormalization scheme [17J, which uses Gp and the physical 
particle masses as basic parameters. Various prescriptions for the renormalization of the 
auxiliary variable tan j3, with specific virtues and flaws, may be found in the literature (for 
a review, see Ref. [H]). For definiteness, we employ the Dabelstein-Chankowski-Pokorski- 
Rosiek (DCPR) scheme [19], which maintains the relation tan/? = V2/V1 in terms of the 
"true" vacua through the condition 5vi/vi = bvijvi, and demands the residue condition 
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Ke'E' AO (M J ^o) = and the vanishing of the A°-Z° mixing on shell as ReS^o^o(M^o) = 0, 
where E^o(g 2 ) and E^o^o^ 2 ) are the renormalized A°-boson self-energy and A°-Z° mixing 
amplitude, respectively. 

The evaluation of the diagrams in Fig. [1] is aggravated by the appearance of three 
7 5 matrices inside closed fermion loops. This leads us to adopt the 't Hooft-Veltman- 
Breitenlohner-Maison (HVBM) scheme [20], which allows for a consistent treatment of 
the Dirac algebra within the framework of dimensional regularization. Then, one has 

75 = ^W7'VW, (2) 
where the totally antisymmetric Levi-Civita tensor is defined in D dimensions as 

!\ if (/i, v, p, a) even permutation of (0,1,2,3), 
-1 if (fi,u,p,a) odd permutation of (0,1,2,3), (3) 
otherwise. 

In fact, we explicitly verified that the naive anticommuting definition of the 75 matrix 
yields ambiguous results, which depend on the way of executing the Dirac traces. Fur- 
thermore, in the renormalization of the pseudoscalar current 

P(x) = Z 2 Z?Z%${x)itfl>{x), (4) 

one needs to introduce a finite renormalization constant Z§, besides the usual fermion 
wave-function and pseudoscalar-current UV renormalization constants Z 2 and Z v of the 
modified minimal-subtraction (MS) scheme, to effectively restore the anticommutativity 
of the 75 matrix [21] . Within QCD, Z\ is known through O(a^) [2l] . Here, we need Zf 
at 0(x t ). We thus need to consider the diagrams depicted in Fig. [2] with the external legs 
amputated, where a cross indicates the insertion of the Fourier transform of P(x) and a dot 
the operator renormalization. Since the A°bb coupling is suppressed in our case, dominant 
contributions only arise from the neutral scalar bosons. Using the mixed commutation 
and anticommutation relations properly distinguishing between 4 and (D — 4) dimensions 
[20] , we decompose the string of gamma matrices appearing in the expression for Fig. Mb) 
into the term proportional to the LO result of Fig. [2](a) that one would obtain with an 
anticommuting 75 matrix and an evanescent remainder, which lives in the unphysical 
(D — 4)-dimensional part of space-time and vanishes in the physical limit D — > 4. Upon 
loop integration, the first term may produce an UV divergence, which would be canceled 
by Z2Z P in Fig. [2](c), while the evanescent remainder may generate an unphysical finite 
contribution to be canceled by Z\ . By explicit evaluation, the latter is found to vanish at 
0(xt), owing to the cancellation of the individual contributions from the z°, h°, H°, and 
A bosons, so that Z\ = 1 for our application. 

We first consider the A —>■ 77 decay. By Lorentz covariance, its transition matrix 
element takes the form 

T = e llvpa e>»{q 1 )er{qz)<f 1 <&A, (5) 
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Figure 2: Feynman diagrams contributing to Zf at Crosses and dots indicate 

the insertions of the Fourier transform of P(x) and its operator renormalization Z2Z P , 
respectively. 



where e^{q) is the polarization four- vector of a photon with four-momentum g M and e^ upcT 
is defined in Eq. (|2D, so that 
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The form factor A is evaluated perturbatively as A = Aq + A as + A Xt + ■ • -. In D 
dimensions, the LO result reads 
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where C = (2 1 ^ 4 /^G^aNc cot (3 with a being Sommerfeld's fine-structure constant, 
Qt = 2/3 is the fractional electric charge of the top quark, 7^ is the Euler-Mascheroni 
constant, and r = M\ /{2m t ) 2 . For r< 1, the function within the square brackets of 
Eq. flZJ) has the expansion 1 + 0(r). The 0(x t ) result A xt = A^J + A xt is composed of 
a counterterm A^J and the contribution A Xt from the proper vertex diagrams in Fig. CD 
We have 
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where 5v/v is the common DCPR counterterm for the two Higgs doublets given in 
Eq. (3.11) of Ref. [22], Ar [17] contains those radiative corrections to the muon lifetime 
which the SM introduces on top of those derived in the QED-improved Fermi model, and 
5m t /m t may be found, e.g., in Eq. (74) of Ref . [23]. In terms of (transverse) self-energies, 
we have 
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We evaluate A® by applying the asymptotic-expansion technique with the help of the 
programs QGRAF [24], q2e, exp [25], and MATAD [26]. Our final result reads 
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77) receives the electroweak 



Comparison with Eq. fl7_D shows that, for r < 1, T(A° 
correction factor [1 — x t (4 + 7/ tan 2 j3)]. 

We now turn to T(A° — > gg). We then need to include the color factor N g /4 = 2 in 
Eq. © and substitute C -> C = {2 l ^/n)G 1 J 2 a s cot (3 and Q t -> 1 in Eqs. © and ©. 
Implementing the appropriate substitutions in the relevant subset of diagrams in Fig. [T] 
and combining the outcome with the counterpart of Eq. (jSJ), the counterpart of Eq. ffTUl) 
is found to be 

a* f-^ - ^1 , (ii) 
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so that T(A° — > (7(7) receives the electroweak correction factor [1 — x t (7 + 10/ tan 2 /?)]. 
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Figure 3: 0(x t ) and 0(a s ) corrections to T(A° — > 77) (a) for M A o = 100 GeV as functions 
of tan/? and (b) for tan/? = 2 as functions of M^o. 

As a check for our computational setup, we also recalculated the 0(a s ) corrections 
to Z\ and, as an expansion in r through 0(t a ), to T(A° — > 77) to find agreement with 
Refs. [21] and [3 El HOl E] , respectively. Notice that the O(r ) term vanishes, so that the 
0(a s ) correction is suppressed for small values of M^o. In fact, as a consequence of the 
Adler-Bardeen theorem [27], the large-m t effective Lagrangian of the A ^^ interaction 
does not receive QCD corrections at any order [2S]. The 0{x t ) and 0(a s ) corrections 
to T(A° — > 77) are compared in Fig. [31 We observe from Fig. [3](a) that the 0(x t ) 
correction amounts to —1.7% at tan/3 = 2 and rapidly reaches its asymptotic value of 
— 1.2% as tan/3 increases, whereas the 0(a s ) correction is positive and independent of 
tan (3. The M^o dependence of the 0(x t ) correction shown in Fig. [3](b) is induced by Ao in 
Eq. (J7J) to which A Xt is normalized. Since it is rather feeble, we may expect the unknown 
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0{r n ) (n = 1, 2, 3, . . .) terms in Eq. (|T0|) to be of moderate size, too. The smallness and 
approximately quadratic M A o dependence of the 0(a s ) correction is due to the absence 
of the leading O(t ) term discussed above. We conclude that the 0(xt) reduction more 
than compensates the 0(a s ) enhancement for M^o <^ 120 GeV. 

The 0(xt) correction to T(A° — > gg) ranges from —2.8% at tan f3 = 2 to the asymptotic 
value —2.1% and partly screens the sizeable 0(a s ) and 0(a 2 s ) corrections of about 68% 
and 23%, respectively, which do have non- vanishing O(t ) terms [8l [T3]. 

In conclusion, we analytically calculated the dominant electroweak two-loop correc- 
tions, of order 0(x t ), to T(A° — > 77), T(A° — > gg), (7(77 — > A ), and cr(gg — > A ) within 
the 2HDM with low- to intermediate-mass Higgs bosons for small to moderate value of 
tan j3 using asymptotic expansion in M\ / (2mt) 2 . To consistently overcome the non-trivial 
75 problem of dimensional regularization, we adopted the HVBM scheme and included the 
finite renormalization constant Z\ of the pseudoscalar current to effectively restore the 
anticommutativity of the 75 matrix. The 0(x t ) term of Z\ was found to vanish. On the 
phenomenological side, the 0(x t ) correction to T(A° — > 77) and (7(77 — > A ) is of major 
importance, since it more than compensates the 0(a s ) enhancement for M^o <^ 120 GeV. 
As for r(v4° — > gg) and u{gg — > A ), the 0(x t ) correction appreciably screens the sizeable 
QCD enhancement, by up to —3%. 

We thank W. Hollik for a useful communication concerning the renormalization of 
tan/3, P. Kant for providing us with the series expansion of Eq. (2.20) in Ref. [TO] , and 
M. Gorbahn, M. Spira, and M. Steinhauser for beneficial discussions. This work was 
supported in part by BMBF Grant No. 05 HT6GUA and DFG Grant No. GRK 602. 



References 

[1] F. Gianotti et al, Eur. Phys. J. C 39, 293 (2005). 

[2] M. Spira, Fortsch. Phys. 46, 203 (1998 ); M.S. Carena and H. E. Haber, Prog. Part. 
Nucl. Phys. 50, 63 (2003); A. Djouadi, [aTXiv:hep-ph/05Mr73"l 

[3] M.M. Muhlleitner, M. Kramer, M. Spira, and P.M. Zerwas, Phys. Lett. B 508, 311 
(2001); D.M. Asner, J.B. Gronberg, and J.F. Gunion, Phys. Rev. D 67, 035009 
(2003); P. Niezurawski, A.F. Zarnecki, and M. Krawczyk, Acta Phys. Polon. B 37, 
1187 (2006). 

[4] P. Kalyniak, R. Bates, and J.N. Ng, Phys. Rev. D 33, 755 (1986); 34, 172 (1986); 
J.F. Gunion, G. Gamberini, and S.F. Novaes, ibid. 38, 3481 (1988). 

[5] J.F. Gunion and H.E. Haber, Nucl. Phys. B278, 449 (1986). 

[6] M. Battaglia and K. Desch, AIP Conf. Proc. 578, 163 (2001); LHC/LC Study Group, 
G. Weiglein et al, Phys. Rept. 426, 47 (2006). 

[7] A. Djouadi, M. Spira, and P.M. Zerwas, Phys. Lett. B 311, 255 (1993). 



8 



[8] M. Spira, A. Djouadi, D. Graudenz, and P.M. Zerwas, Nucl. Phys. B453, 17 (1995). 

[9] M. Spira, A. Djouadi, D. Graudenz, and P.M. Zerwas, Phys. Lett. B 318, 347 (1993). 

[10] R. Harlander and P. Kant, JHEP 0512, 015 (2005). 

[11] U. Aglietti, R. Bonciani, G. Degrassi, and A. Vicini, JHEP 0701, 021 (2007). 

[12] V. Ravindran, J. Smith, and W.L. van Neerven, Nucl. Phys. B704, 332 (2005). 

[13] K.G. Chetyrkin, B.A. Kniehl, M. Steinhauser, and WA. Bardeen, Nucl. Phys. B535, 
3 (1998). 

[14] R.V. Harlander and W.B. Kilgore, JHEP 0210, 017 (2002). 
[15] R.V. Harlander and F. Hofmann, JHEP 0603, 050 (2006). 
[16] J. Brod, F. Fugel, and B.A. Kniehl, in preparation. 

[17] A. Sirlin, Phys. Rev. D 22, 971 (1980); W. Hollik, E. Kraus, M. Roth, C. Rupp, 
K. Sibold, and D. Stockinger, Nucl. Phys. B 639, 3 (2002). 

[18] A. Freitas and D. Stockinger, Phys. Rev. D 66, 095014 (2002). 

[19] PH. Chankowski, S. Pokorski, and J. Rosiek, Nucl. Phys. B423, 437 (1994); 
A. Dabelstein, Z. Phys. C 67, 495 (1995). 

[20] G. 't Hooft and M.J.G. Veltman, Nucl. Phys. B44, 189 (1972); P. Breitenlohner and 
D. Maison, Commun. Math. Phys. 52, 11 (1977). 

[21] S.A. Larin, Phys. Lett. B 303, 113 (1993). 

[22] A. Dabelstein, Nucl. Phys. B456, 25 (1995). 

[23] M. Butenschdn, F. Fugel, and B.A. Kniehl, Nucl. Phys. B772, 25 (2007). 
[24] P. Nogueira, J. Comput. Phys. 105, 279 (1993). 

[25] R. Harlander, T. Seide nsticker, an d M. Steinhauser, Phys. Lett. B 426, 125 (1998); 
T. Seidensticker, |hep-ph/9905'298l 

[26] M. Steinhauser, Comput. Phys. Commun. 134, 335 (2001). 

[27] S.L. Adler and W.A. Bardeen, Phys. Rev. 182, 1517 (1969). 

[28] B.A. Kniehl and M. Spira, Z. Phys. C 69, 77 (1995). 



9 



